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ABSTRACT 

The local cosmic-ray (CR) spectra are calculated for typical characteristic regions of a cold dense 
molecular cloud, to investigate two so far neglected mechanisms of dust charging: collection of 
suprathermal CR electrons and protons by grains, and photoelectric emission from grains due to the 
UV radiation generated by CRs. The two mechanisms add to the conventional charging by ambient 
plasma, produced in the cloud by CRs. We show that the CR-induced photoemission can dramatically 
modify the charge distribution function for submicron grains. We demonstrate the importance of the 
obtained results for dust coagulation: While the charging by ambient plasma alone leads to a strong 
Coulomb repulsion between grains and inhibits their further coagulation, the combination with the 
photoemission provides optimum conditions for the growth of large dust aggregates in a certain region 
of the cloud, corresponding to the densities n(H 2 ) between ~ 10 4 cm -3 and ~ 10 6 cm" 3 . The charging 
effect of CR is of generic nature, and therefore is expected to operate not only in dense molecular 
clouds but also in the upper layers and the outer parts of protoplanetary discs. 

Subject headings: ISM: dust ISM: clouds - ISM: cosmic rays 


1. INTRODUCTION 


Interstellar dust grains in dense molecular clouds 
are s ubject to several elect ri c charging pr o cesses 
fe.g.. IDraine fc Salneted 119791 IDraine fc Sutinl 119871 : 
iWeineartner fe Draindl2001bfl . The resulting net electric 
charge carried by micron or sub-micron size grains has 
important consequences for the chemical and dynamical 
evolution of mole cular clouds: it affects the process o f 
dust coagulation (|Okuzumil 120091 IDominik et al.l [2007h . 
the rate of grain-catalyzed electron-io n recombination 
(|Mestel fc Spitzeil 119561 : iWatsonl I1974D. the amount of 
gas-phase elemental depletion dSoitzeil 11941 ), and the 
electrical resistivity of t he cloud’s plasma ( Elmegreenl 
119791: IWardle fc Ng)ll999fl . The resistivity, in turn, con¬ 
trols the coupling between the neutral gas and the in¬ 
terstellar magnetic field, and eventually the dynamics of 
gravitational colla pse of molecular clouds and the form a- 
tion of stars fe.g.. iNakano et aLl '2002: Sh u et all 1200611 . 


Collisions of dust grains with the plasma of thermal 
electrons and ions from the gas (hereafter, cold plasma 
charging) represent an im portant dust charging process 
in m olecular clouds fe.g.. IDraine fc Sutinl 119871 : iDrainel 
l2onl) . Since electrons of mass m e have a thermal speed 
which is much larger than that of ions of mass nrii (by 
the factor \Jmi/m e 1), grains acquire by this pro¬ 
cess a (predominantly) negative charge. The photo¬ 
electric effect (also called photoemission), on the other 
hand, results in positive charging of dust grains, and 
is set by the radiation field in the cloud at energies 
above a few eV. Photoemission is an important charg¬ 
ing proce ss for diffuse gas with visual extinction Ay < 
10_Je.g., iBakes fc Tielensl 11994 iWeingartner fc Drainel 
l2001blj . As the interstellar radiation field is exponen¬ 
tially attenuated with increasing Ay, photoemission is 
usually neglected to compute the charge distribution of 
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grains in the dense ga.s of molecular cloud c ores (e.g., 
lUmebavashi fc Nakanol[l980t iNishi et ahll 19911) . 

Cold plasma charging and photoemission are usually 
assumed to be the dominant grain charging mechanisms 
in the cold interstellar medium. In this paper we study 
the effects of cosmic rays (CRs) on the charging of sub¬ 
micron dust grains in molecular clouds. We focus on 
two charging processes that contribute in addition to the 
cold-plasma charging, but have been neglected so far. By 
calculating the local CR spectra for typical cloud regions, 
we investigate the effects of (?) collection o f suprat hermal 
CR electrons and protons by grains dShchekinovt 1200 71 
and (u) photoelectric emission from grains due to the 
UV radiation generated by CRs. Using the cold-plasma 
collection as the “reference case”, we show that the pho¬ 
toelectric emission can dramatically modify the charge 
distribution function for dust in almost the entire cloud, 
and discuss important implications of the obtained re¬ 
sults. In particular, we point out that while the cold- 
plasma charging alone leads to a strong Coulomb repul¬ 
sion between grains and inhibits their further coagula¬ 
tion, the combination with the CR-induced photoemis¬ 
sion provides optimum conditions for the growth of large 
dust aggregates in a certain region of the cloud. 

2. CR PROPERTIES RELEVANT TO DUST CHARGING 

The specific intensities (or spectra) of CR protons and 
electrons inside a dense molecular cloud are determined 
by the interstellar CR spectra. In order to constrain 
the trend of the interstellar spectra at high energies 
{E > 500 MeV), we use the latest results of the Al¬ 
pha Magnetic Spectrometer (AMS-02), mounted on th e 
International Space Station (jAguilar et a,l .1 IaTmI 12015(1 . 
The high-energy spectrum slope is —3.2 for electrons, 
while for protons it is —2.7. 

At lower energies, the shape of the interstellar CR spec¬ 
trum is highly uncertain due to the effect of Solar mod- 
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ulation (see, e.g.. iPutze et al.lf201lD . Data collected by 
the Voyager 1 spacecraft from a region beyond the Solar 
terminat ion shock are extremely useful in this context 
fWebheil 119981) . as they provide a lower limit to the CR 
spectrum in the energy range from ~ 5 MeV to ~ 50 MeV 
that should be not too different from the actual inter¬ 
stellar value dStone et all 1201 ill ). However, two caveats 
should be kept in mind. First, Voyager 1 has not yet 
entered interstellar space, as the magnetic field detected 
by magnetometers on board the spacecraft still retains 
some characteristics o f the Solar wind magnetic field 
(jBurlaga fc Nessi [2014t h Therefore, the measured pro¬ 
ton fluxes may still contain a fracti on of anomalous CRs 
from inside the heliosphere dScherer et alJl2008h . Second, 
even if Voyager 1 were in the interstellar space, there is no 
guarantee that the measured spectra of CR protons and 
electrons are representative of the average Galactic spec¬ 
tra, because the contribution of local sources (i.e. within 
parsecs from the Sun) is difficult to quantify. Neverthe¬ 
less, the Voyager 1 data provide the only direct obser¬ 
vational constraint presently available on the low-energy 
spectra of CRs, and therefore they cannot be discarded. 
In particular, we use the results of [Stone et all (2013), 
obtained from data collected by Voyager 1 since August 
2012, when the spacecraft was at an heliocentric distance 
of 122 AU. 

We model the low-energy behavior of the proton and 
electron spectra with power-law dependencies. To de¬ 
scribe a crossover to the high-energy scalings, we employ 
the simple analytical expression 

jk{E) = C (E + E 0 )P ^cm-V'sr- 1 , (1) 

where k = p,e and the crossover ener gy Eg = 500 MeV i s 
the same for both species. Following iStone et all (1 201 3 l'l. 
for CR electrons we adopt a low-energy spectrum slope 
of a = —1.5, while for protons we explore two extreme 
cases with a = —0.8 and a = 0.1. The resulting inter¬ 
stellar spectra are presented in Fig. [T] other parameters 
for Eq. (|T]) are listed in Table [1] Combinations of the two 
proton spectra with the electron spectrum are termed be¬ 
low as the model Jf’ (“High”, a = —0.8) and model Jz? 
(“Low”, a = 0.1). 


TABLE 1 

Parameters of interstellar CR spectra, Eq. {TJ. 


species ( k ) 

c 

a 

P 

(3 — a 

electrons 

2.1 x 10 18 

-1.5 

1.7 

3.2 

protons (model Jtf?) 

2.4 x 10 15 

-0.8 

1.9 

2.7 

protons (model 

2.4 x 10 15 

0.1 

2.8 

2.7 


The number densities for the CR species are given 
by n C R,k = 4tt f™ cut dE j k (E)/v k (E), where v k (E) is 
the velocity for the energy E. In order to preserve 
charge neutrality of CRs, we set ncR.p = ncR.e; which 
yields the lower energy cutoff E cut — 5 keV and ~ 
2 MeV for the models J4? and „Sf, respectively. The 
corresponding energy densities are defined by £cn,k = 
An^dE Ej k (E)/v k (E). The total energy density, 



Fig. 1.— Interstellar spectra of CR protons (a) and electrons 
(b). The black solid lines show the parameterized s pectra given 
by Eq. Q. open circles represent the Voyager 1 data tStone_et_aL| 
2012), gray dashed lines - the AMS-02 data 1 A"i]i 1 ar et al .11 201 41 . 
2015 1. A combination of the given electron spectrum with the 
proton spectrum for a = —0.8 and a = 0.1 is referred to as the CR 
model Jif and JC, respectively. The respective lower energy cutoff 
E cut (— 5 keV and ~ 2 MeV) is indicated by the gray vertical bars. 

£cr,p + £CR,e, is dominated by protons and varies be¬ 
tween ~ 0.78 eV cm -3 (model Jzf) and ~ 1.54 eV cm -3 
(model Jif). 

We note that our choice of a single interstellar spec¬ 
trum for CR electrons and two possible spectra for CR 
protons is arbitrary. As discussed above, we adopt the 
view that the Voyager 1 data represent lower limits to 
the actual spectra, due to residual modulation of the in¬ 
terstellar fluxes at the current position of the spacecraft. 
For simplicity, this remaining uncertainty is attributed 
to the protons only: We select the “minimum” proton 
spectrum (a = 0.1) compatible with the data, and the 
“maximum” spectrum (a = —0.8) providing the upper 
bound for the available observational data on the CR 
ionization rate (see Sec. EH- 

2.1. Local CR spectra 

In order to concentrate on the charging effects induced 
by CRs, we adopt an idealized ID (slab) model of a 
den se core em bedded in a molecular cloud (for details, 
see lPadovani et aD 120091 ). We assume that CRs propa¬ 
gate normally to the surface in straight lines, with a half 
of the interstellar CR flux incident on each side of the 




































3 


O X c 



Fig. 2.— Local CR spectra for the model (solid lines) and model jSf (dashed lines). The three panels show the results for the three 
typical characteristic regions of a dense core: “outer” ( O ), “inner” ( X ), and “center” (C). The thin lines represent a cold Maxwellian plasma 
background, the thick lines are for the suprathermal propagated CR spectra. The gray and black lines correspond to electrons and protons, 
respectively. The presented results are for the case when H + are the dominant ions in a cold plasma; when heavier ions dominate, the 
Maxwellian spectra for ions should be divided by a square root of the corresponding atomic mass number. 



Fig. 3.— (a) CR ionization rate ( as a function of the column 
density of molecular hydrogen 7V(H2) for the model Jj? (black line) 
and model (gray line), the densities corresponding to the regions 
O, I , an d C are indicated. The sy mbols are the obse rvational val¬ 
ues fr om Indriolo & McCall (201 2), gray solid cir c les) ■ ICaselli et al.l 
1119981 . black solid circles), and IMaret fe Berginl (120071 . empty cir- 
cle). (b) Comparison of the attenuated interstellar flux of UV pho¬ 
tons (thick lines) with the local UV flux generated by CRs (thin 
lines). The thin black and gray lines represent the CR-generated 
fluxes for the model Jif and .2”, respectively. The dashed and solid 
lines show the results for Ry = 3.1 and 5.5, respectively. 


cloud. 1 This simple model, neglecting integration over 
the incidence angles, yields the local spectra that ap¬ 
proach the exact results in the inner core region, while in 
the outer region the fluxes are only slightly overestimated 
(by less than 30%, details will be presented elsewhere). 

For the core we use the d e nsity profile of a Bonnor- 
Ebert sphere dBonnoil 119561 : lEberti I1955H . which typ- 
icall y well reproduces observations of starless cores 

(e.gJAlves et alJ200H iKeto fc Casellill2008l: lAndre et al.1 

I2014D . In particular, we consider a centrally concentrated 
core, such as L1544, where the volume density within 
the central 500 AU is n( H 2 ) = 2 x 10 7 cm~ 3 , one of th e 
largest known for starless cores (|Keto fe Casellil I2010D . 
The molecular core has a radius of about 0.1 pc, be¬ 
yond which the density drops below n( H 2 ) = 1 0 4 cm - 3 
and photoprocesses becomes important (IKeto fe Casellil 
12010 ). We further assume that the core is embedded 
in an envelope of lower density, to simulate the location 
of the dense core in a molecular cloud such as Taurus. 
Such envelope extends up to a distance of 4.4 pc from the 
core center, and the average density within a radius of 
4.4 pc is 230 cm -3 . With these properties, the cloud has 
a mass column densit y of 144 pc~ 2 , which is typical 

of molecular clouds dRoman-Duval et ahll2010D . 

To demonstrate the extent to which the dust charges 
in dense molecular clouds are expected to be affected by 
CRs, we consider three characteristic regions of the em¬ 
bedded core: the outer boundary [n(H 2 ) = 10 4 cm -3 , 
7V(H 2 ) = 3.2 x 10 21 cm -2 ], the inner core [n(H 2 ) = 
10 6 cm -3 , iV(H 2 ) = 2.8 x 10 22 cm -2 ], and also the core 
center [n(H 2 ) = 2 x 10 7 cm -3 , JV(H 2 ) = 10 23 cm -2 ]. We 
shall refer to these regions as “outer” (O), “inner” (X), 
and “center” (C), respectively. 

The CR protons and electrons penetrating the cloud 
ionize the neutral gas, i.e., produce a local plasma envi- 


1 More advanced models should take into account the fact that 
molecular clouds are magnetized and CRs gyrate along magnetic 
field lines, in addition of being scattered by magnetic fluctuations 
on the scale of t he particle gyrorad i us. F o r detailed treatment of 
these effect s, see IPadovani &: Gall: (1201 11 1. IPadovani et al.1 (120131 . 
1201411 , and IMorlino7r*^ab!cil ~ il20TE ). 
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ronment whose properties are determined by the value 
of iV(H 2 ). On the other hand, the ionization yields a 
major contribution to the energy loss of CRs (see Ap¬ 
pendix m, thus modifying_the interstellar spectra to 
the local form dPadovani et al.l[2009il . The low-energy 
CR species as well as the multiple generations of elec¬ 
trons and ions formed in the ionization avalanche (see 
Appendix 03 rapidly lose their energy: the character¬ 
istic energy loss timescale T st0 p, given by Eq. (IB2I) . is 
always much shorter than the timescale of recombina¬ 
tion. Therefore, the steady-state energy distribution of 
electrons and ions in a dense cloud can be viewed as a su¬ 
perposition of two distinct parts: a cold Maxwellian peak 
where electrons and ions accumulate and eventually re¬ 
combine - that dominates the total plasma density, and a 
suprathermal tail representing the modified (propagated) 
CR spectrum - that determines the ionization rate. 

The local steady-state equilibrium is determined 
by the balance between the CR ionization of 
H 2 an d var ious reco mbination p rocesses (see, e.g., 
iQppenheimer fe Dalgarnol 11974 lMcKedll989lh Acorn- 
petition between the dissociative recombination with 
molecular ions and the radiative recombination with 
heavy metal ions (in the presence of charge-transfer re¬ 
actions), occurring in different core regions, can sig¬ 
nificantly alter the magnitude of the electron fraction 
x e = n e /n{ H 2 ), in particular modify the dependence on 
rc(H 2 ). For the L15 44 core, we e mploy the following in¬ 
terpolation formula dCaselli et al.ll2002 . their model 3): 


x e ~ 6.7 x 1CT 6 




( 2 ) 


where £ is the CR ionization rate. 2 In the following, we 
assume that the electron density n e is equal to the den¬ 
sity of all ion species formed in the ionization avalanche 
(i.e., the contribution of charged grains into the charge 
neutrality is negligible, see Sec. 

Figure [2] shows the combination of the propagated 
spectra and the cold Maxwellian background. The en¬ 
ergy E[nt at which the intersection between the two 
curves occurs is practically the same for all considered 
cases, Flint — 1.5 x 10“ 2 eV. 


TABLE 2 

CR IONIZATION RATE £ (S -1 ) AND PLASMA ELECTRON FRACTION 
X e , FOR THREE CHARACTERISTIC CORE REGIONS (O, X, AND C) 
AND TWO CR MODELS (_Sf AND . 



c 


X e 



model 

model 

model 

model 

o 

6.0 x 10~ 1Y 

6.5 x 10~ ie 

9.4 x 10" 8 

3.1 x 10~ 7 

X 

5.6 x 10- 17 

3.0 x 10~ 16 

6.9 x 10" 9 

1.6 x 10 -8 

c 

5.3 x 10- 17 

1.8 x lO" 16 

1.3 x 10 -9 

2.3 x 10" 9 


Figure [3ji shows the ionization rate as a function of 
the molecular H 2 column density for the two CR models. 
The characteristic values of £ and the corresponding x e 
are summarized in Table [2j For comparison, the values 
of ( obtained by integrating the Voyager 1 fluxes down 

2 We note that the results presented in Sec. [3] practically do not 
depend on the precise form of the formula for x e . 


to the lowest measured energy, without any extrapola¬ 
tion, are £ = 1.2 x 10 -17 s -1 and £ = 2.8 x 10 -18 s _1 
for protons and electrons, respectively. These should be 
considered as lower limits to the interstellar value of £. 
Notice that the contribution of CR electrons to the av¬ 
erage interstellar CR ionization rate, often neglected in 
the past, could be significant (although in general it is 
smaller than that of protons and heavier nuclei). 

We point out that the main positive charge carrier 
changes across the cloud, because the molecular freeze- 
out becomes more effici ent toward the core center (see, 
e.g., iTafalla et al. 2002). In this paper we do not dis¬ 
cuss a complex plasma chemistry, but instead consider 
two extreme cases when the dominant ions are either H + 
(for the strongly depleted inner dense regions) or HCO + 
(for the outer regions). In Sec. 13.31 we demonstrate that 
the plasma composition has only a minor effect on the 
obtained results. 

2.2. Local radiation field 

The local radiation is generated by CRs via the fol¬ 
lowing three (prime) mechanisms: The CR electron 
bremsstrahlung, the 7r° decay, and the H 2 fluorescence. 

We ca lcu late t he br emsst rahlung spectrum following 
iBlumenthal fc Gouldl (1970). For all three characteris¬ 
tic core regions, the resulting photoemission flux from 
a grain is much smaller than the collection flux of the 
surrounding cold-plasma ions, which indicates that the 
CR brem sstr ahlung cannot contribute to dust charging 
(see Sec. 13.21) . For the photon_spectrum due to 7r° de¬ 
cay, we follow iKamae et ahl ( 2006 b This yields the pho¬ 
toemission flux which is even smaller than that due to 
bremsstrahlung and, hence, is negligible, too. Note that 
for energies larger than the pion production threshold 
(280 MeV), both our CR proton spectra coincide and 
remain unmodified up to 7V(H 2 ) ~ 10 23 cm - 2 . 

Fina lly, following Eq. (21) in lCecchi-Pestellini fe Aiellol 
GMD, we compute the H 2 fluorescence generated by 
CRs in the Lyman and Werner bands. The resulting 
flux of UV photons Fuv (in the energy range between 
11.2 and 13.6 eV) can be approximately calculated as¬ 
suming that the band excitation rates, being normalized 
by £, do not depend on the shape of the CR spectrum: 3 

Fdv=!96 „(_1_)(_C_) (3) 

/ AWMv WfivV 5 cm - v , 
£10 21 cm 2 mag 1 ) \3.2/ 

Here, w is the dust albedo at ultraviolet wavelengths, and 
Ry is a measure o f the slope of the extinction at visible 
wavelengths (e.g., lDraineli201l|). Ass uming lj = 0.5 and 
Ry = 3.1 (|Cecchi-Pestellini fe Aiellolll992f) . 4 and substi¬ 
tuting the typical gas-to-extinction ratio of V(H 2 )/Ay = 

3 In fact, the excitation and ionization cross sections have some¬ 
what different dependencies on energy. For this reason, the nor¬ 
malized excitation rate by CR electrons, obtained with the spectra 
sh own i n Fig. m is about 2.2-2.8 t imes larger than that reported 
in lCecchi-Pestellini fc Aiellol (| 19921 ). However, we were not able to 
find a reliable expression for the excitation cross section by protons, 
and therefore employ the approximate Eq. 

4 According to I Gordon! (120041) , in dense cores with Ry > 5 the 
dust albedo may be closer to 0.3. 
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10 21 cm 2 mag 1 we get Fuv — 1830(^/10 17 s 7 ). In 
Sec. [3] it is shown that the resulting photoelectric emis¬ 
sion from dust grains can significantly exceed the cold-ion 
collection. Thus, the strongest radiation field generated 
by CRs is due to the H 2 fluorescence. 

The interstellar radiation field is exponentially attenu¬ 
ated in the cloud. The specific intensity at the frequency 
v decreases as I(is) = I is (is)e ~ Tu , where Iis(is) is the 
i ntensit y of the interstellar radiation field according to 
iDrainel (|2011f ) and t v = H„/1.086 is the optical depth. 
By integrating I (is)/his over his > 10 eV (at lower fre¬ 
quencies, the photoelectric yield from dust grains rapidly 
falls off, so the radiation does not contribute to the pho¬ 
toemission), we derive the interstellar photon flux as a 
function of iV(H 2 ). The latter is plotted in Fig. [3jr and 
compared with F\jy vs. iV(H 2 ) calculated for the models 
Jif and jSf. We see that the photon flux generated by 
CRs only slightly decreases (model sZC) or remains prac¬ 
tically constant (model FZ) in the shown range of N (H 2 ); 
even for the maximum value of Ry = 5.5, the interstellar 
flux becomes negligible at 1V(H2) > 3 x 10 21 cm -2 . We 
conclude that for all three core regions the radiation field 
is solely due to CRs. 

3. CHARGE DISTRIBUTION FUNCTION OF DUST 

The discrete charge distribution N(Z) = Nz for dust 
grains of a given size a is normalized to the total differen¬ 
tial dust density at that size, i.e., Yz^z = drid(a) / da. 
The charge distribution is derived from the detail ed e qui¬ 
libriu m of the charging master equation dDraine fe Sutinl 

sm 


Je(Z + l)Nz+l 


J2MZ) + J PE (Z) 


Nz- 


( 4 ) 


Here, J e ,i(Z) = J^(Z) + J^f-(Z) is the electron/ion 
collection flux, which has the contributions from cold 
Maxwellian plasma background (first term) and from 
suprathermal low-energy part of the CR spectra (second 
term), and J PE (Z) is the photoemission flux due to the 
local radiation field. 

We introduce the floating potential ipz = Ze 2 /a of a 
particle with the charge Ze. The collection fluxes of cold 
electrons and ions are obtained from the so-calle d “or¬ 
bital motion limited” ( OML) approximation (WhioDlc 
fTMl IFortov~i r ani2005h : 


J^(Z) = 2\[Yk a 2 n e v e 
and 

J™(Z) = 2V2na 2 mvi 


gV)z/feBT 

Z < 0; 

(1 + p z /k B T), 

z> 0 , 

(1 - p z /k B T), 

Z< 0; 

g—Vz /ksT 

Z>0, 


( 5 ) 


( 6 ) 


where u e y = \Jk B T / m e> i are the thermal velocity scales. 
Note that for the sake of clarity the (minor) effect of the 
polarization interaction, omitted in these equations, is 
discussed later. 

The relation between the electron and ion densities is 
determined from the charge neutrality condition, 

n e = '^rii + ( Z)n d , 


where (Z)n d = / a amax da Yz ZNz(a) is the charge num¬ 
ber density carried by dust. In starless dense molecu¬ 
lar clouds studied in this paper, the dust contribution is 
negligible as long as the ratio n(H 2 )/C is below a cer¬ 
tain value - e.g., for single-size grains of a = 0.1 pm 
t he neutr al densi t y sh ould be less than ~ 10 11 cm~ 3 
(lUmebavashi fe Nakanc I1990H . for the MRN distribu¬ 
tion ( Mathis e t_ah|11977) the ap proach works at n(H 2 ) < 
10 7 cm -3 (|Nakano et al.ll2002f) . Therefore, for the cal¬ 
culations below we set n e = Yi n i f° r all three cloud 
regions, which implies that our results are independent 
of details of the grain size distribution and demonstrate 
the generic CR effects on dust charging. 

We note that in situations where the dust contribu¬ 
tion to the charge neutrality is not negligible, e.g., at 
higher n( H 2 ), the distribution of grain charges is de¬ 
termined by the particular form of dn d (a)/da which is 
quite uncertain inside dense molecular clouds (see, e.g., 
IWeingartner fe Drainell2b01atlKim et al.lfl994T ). Some in¬ 
sights into possible forms of the size distribution in such 
environm ents can be gain ed from numerical simulations 
(see, e.g.. lOrmel et alJl2009H . The role of dust becomes 
particularly important in dense cores containing proto- 
stellar sources, where small grains (the main carriers of 
negative cha rge) a re a bundantly produced due to shock 
shattering (IGuillet et al.ll201lh . 


3.1. “Reference case”: Cold-plasma charging 

In the usual approach, the dust charging in cold molec¬ 
ular clouds is caused by collection of electrons and ions 
from a cold plasma background. The charge distribution 
in this case is readily obtained by substituting Eqs. m 
and © in Eq. 


iV+i 1 

No Vm( 1 + ip) ’ 
iV_ 1 \/ni. 

No 1 + ’ 

N —2 

“ 1 + 2Cp ’ 


( 7 ) 


The distribution depends on two dimensionless numbers: 
The normalized floating potential of the unit charge (or 
the inverse normalized temperature) (p = e 2 /ak^T, and 
the effective ion-to-electron mass ratio to (or the effec¬ 
tive atomic mass number A) determined by the partial 
contributions of all ions: 

_ _ 1.67 42.8 _ 1 ^sr^rii 1 

^ a^Tio’ \ffli \[A ^— J n e \J~Ai 

where a^ is in units of pm, T\o is in units of 10 K, and 
Ai is the atomic mass number of the zth ion species. 
Both numbers are large: (p ~ 17 for T = 10 K and 
the largest grains of the MRN distribution (a ~ 0.1 pm), 
while to = 1836 for a hydrogen plasma. This implies that 
(i) the abundance of the Z = —2 state is exponentially 
small, and (ii) one can neglect unity in the denomina¬ 
tors of Eq. ©, so the ratios N±i/N 0 oc aT have a simple 
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universal dependence on temperature and dust size. Fur¬ 
thermore, since N-i/N+i = to, the negatively charged 
state is at least three orders of magnitude more abundant 
than the positive one. 

Thus, submicron grains in a cold plasma are either neu¬ 
tral or singly negatively charged. For a given tempera¬ 
ture, most of the dust smaller than a^ ~ 1.67/(\/mXio) is 
neutral, while larger grains are mostly negatively charged 
(for a hydrogen plasma at T = 10 K, the transition oc¬ 
curs at a ~ 400 A). This simple scaling holds as long as 
(p 1, i.e., for < Tf (j . Larger grains/aggregates (in 
the micron-size range, or if the temperature increases) 
become multiply negatively charged; the charge distri¬ 
bution remains narrow, with the average charge ~{Z) 
about a few (p~ x . 


3.2. Effect of CRs 

Figure [2] demonstrates that the CR proton spectrum 
at lower energies is by many orders of magnitude lower 
than the electron spectrum. Therefore, the former gives 
a negligibly small contribution to the collection flux, as 
compared to the flux of CR el e ctron s. Following the ap¬ 
proach bv IDraine &; Salneterl (|1979D . we calculates the 
electron collection flux as 

pOO 

J e CR ~ ttu 2 / dE 4 tt j e (E) [s e (E ) - 6 e (E)\ , (8) 

J Eint 

where s e (E ) is the sticking coefficient (probability) for 
impinging CR electrons and S e (E ) is the yield of the sec¬ 
ondary electrons (see details in Appendix [C]). The lower 
limit of integration is equal to the intersection energy 
E- lnt ~ 1.5 x 10 -2 eV, at which the CR electron spectra 
in Fig. [2] cross the Maxwellian curves. 

Note that in Eq. ©, under the integral we o mitted 
the OML factor (1 + tpz/E) (see, e.g., Horanvi et al.1 
I1988H . which determines the charge dependence of the 
flux (and leads to Eq. (© for the Maxwellian spectrum). 
As shown below, the charge distribution is usually con¬ 
centrated within \Z\ < 3; since the magnitude of the 
floating potential is e 2 /a ~ 10 meV, while the energy of 
CR electrons contributing to the charging is E ~ 10 eV, 
the OML correction can be safely neglected and hence 
J,F R does not depend on Z. 

The photoelectric emission flux from a grain J FE (Z) is 
determined by the total specific intensity I{v) of the local 
radiat ion field, and depends on the photoemission yield 
Y(u) (W eingartner et aD 120061) and t he absorption effi¬ 
ciency QabslA) ( Draine fe Had 1200% iDra.inel 120 1 ll) . As 
shown in Fig. ©>, the interstellar UV field (integrated 
over the frequency range of 10 eV <hv< 13.6 eV) is ex¬ 
ponentially attenuated with iV(H 2 ), and so for all three 
core regions the local radiation is solely due to the H 2 
fluorescence generated by CRs. The resulting photoe¬ 
mission flux can be approximately calculated as 

J PE = ira 2 J dv ^ )Q aM 

— KO, 2 F V y (Y(v)Q ahs (v)) V v, 

where F uv is given by Eq. © and (F(i/)Qabs(u))uv is 
averaged over the Lyman and Werne r bands. The l atter 
value may vary for different materials (IWeingartner et al.l 
2006), but this variation is not very significant; for the 


calculations below we set (E(^)Qabs(i / ))uv = 0.2, which 
corresponds to a carbonaceous grain of a ~ 0.1 fim. Here, 
the dependence of J PE on Z is neglected too, since photo¬ 
electrons have a broad energy distribution (roughly lim¬ 
ited from above by the energy of the UV p hotons minus 
the work function, see, e.g.. iDrainel 11978( 1 and, hence, 
the floating potential cannot noticeably affect the pho¬ 
toemission flux. 

Thus, the CR effect on the charge distribution is de¬ 
termined by the collection flux of CR electrons J ER and 
the photoemission flux J PE . These fluxes compete, re¬ 
spectively, with the electron and ion fluxes of a cold 
plasma background and, hence, their relative magnitude 
can be conveniently quantified by the dimensionless num¬ 
bers ccr and e PE : 

J e CR = e C R J e M (0), J PE = ep E ^ 0). 


Both numbers do not depend explicitly on a [a weak im¬ 
plicit dependence on dust size is via s e (E , a) and S e (E, a) 
for ccr and via Y(is, a) for e PE ]. From Eqs. © and © it 
follows that J e M (0) oc n e VT and ^VJ^O) oc n e ^T/A, 
and so the relations to the plasma parameters and the 
UV flux are given by 


ecR oc 


1 

n e y/T ’ 



(9) 


The abundances of the singly charged states are 
straightforwardly derived from Eq. 


N + i _ 1 + e PE 

No y/l%(fp + CCR) ’ 

N- 1 _ \Z?fi( 1 + ccr) 
No ip + e PE ’ 


( 10 ) 

( 11 ) 


for the multiply charged states the following recurrent 
relations are obtained: 


Z >2: 


-Z> 2 : 


N+z _ e PE 
N +(z- 1 ) \/fti(Z(p + ccr) ’ 
N -\z\ _ \AfiepR 

N_(|z|_i) \Z\(p + e PE 


( 12 ) 

(13) 


Following the estimates in Sec. 13.11 the unity was ne¬ 
glected with respect to (p in these equations (as one can 
see by comparing with Eq. 0) ; for the multiply charged 
states, we also omitted the exponentially small terms, 
assuming that ccr, e PE e~ v . 

The charge distribution is slightly modified when the 
polarization interactions are taken into account (see, 
e.g., Eqs. (3.3) and (3.4) in IDraine fe Sutiiil 119871) : In 
Eqs. (TT01) and (fill) , one has to add the term ybn^/2 
to the unity in the numerator, and multiply (p in the 
denominator with 2. For the multiply charged states 
(both positive and negative), one needs to replace \Z\ in 
the denominator of Eqs. m and (THU) with the product 

\/W\(\fW\ + !)■ 

One could also take into account the fact that CR elec¬ 
trons (contribution to the collection flux J^ R ) have some 
finite effective temperature Xcr. In this case, the collec¬ 
tion flux for Z < 0 and, hence, the r.h.s. of Eq. m has 
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Fig. 4.— Effect of CRs on the dust charge distribution in dense molecular clouds. The distribution function Nz is described by Eqs. (1101 )- 
im the results for (f) = 17 and A = 29 (HCO - * - ions) are plotted for different CR models (Jf and Jt?) and different characteristic regions 
(O, X, and C) of the core. The left panel shows the results for typical parameters determining the magnitude of epE, in the right panel 
epE is increased by a factor of 10. The insets in the lower row depict the “reference case”, described by Eq. 0, where the CR effects are 
neglected. The uncharged state Z = 0 (shaded) is normalized to unity. 


to be multiplied with the OML factor e Z I^ CR , w l 16 re 
<£cr = e 2 /ak B T C R = f{T/T c r). As we pointed out 
above, the energy of CR electrons contributing to the 
charging is of the order of 10 eV, so <^cr is very small 
(< 10 -2 for the parameters used above) and the resulting 
effect is indeed negligible. 

TABLE 3 

Dimensionless numbers ccr and £ pe characterizing the CR 
EFFECT. The values correspond to conditions for the left 
panel of Fig. [4] For the right panel, e PE should be 
MULTIPLIED BY 10; FOR DIFFERENT IONS IT SHOULD BE RESCALED 
AS epE OC VA. 


eCR _ _ tpE 



model ££ 

model 

model £? 

model 

o 

2.67 x 10“ 3 

9.33 x 10 -3 

2.38 X 10 2 

8.71 x 10 2 

X 

3.35 x ltr 4 

1.34 x 10“ 4 

29.7 

74.3 

c 

8.56 x 10“ 5 

4.31 x 10“ 5 

7.55 

15.0 


Figure U illustrates the CR effect on the charge distri¬ 
bution in different regions of the cloud. The left panel 
shows Nz calculated for typical plasma density n e [from 
Eq. ©] and parameters determining the magnitude of 
the CR-generated UV flux F B y [Eq. d3j with cu = 0.5 
and Ry = 3.1]. However, both n e and F\jy are known 
only approximately: The ionization fraction might be a 
factor of 2-3 lower t han that give n by Eq. ([2]) (see, for 
example, model 1 of iCaselli et al.ll2002fl . while the UV 


flux is approximately doubled for Ry > 5 (typical for 
dark dense clouds, see also Footnote SI). Furthermore, 
as pointed out in Footnote [3l Fuv can exceed the value 
given by Eq. f3j due to the difference between the ioniza¬ 
tion and excitation cross sections. As one can see from 
Eq. ©, a combination of all these uncertainly factors 
may easily increase the value of epE by an order of mag¬ 
nitude. Therefore, in the right panel we show Nz for 
epE multiplied by 10 (and otherwise the same parame¬ 
ters). The role of CRs in dust charging becomes partic¬ 
ularly evident when the obtained results are compared 
with the “reference” charge distribution depicted in the 
insets, where the CR effects are neglected. 

From Eqs. m and (fl3l) we infer that, generally, the 
asymptotic form of Nz at large \Z\ is a Poisson dis¬ 
tribution P(Z\ A): For positive charges, it operates at 
Z > £cr/) 5 and is determined by the Poisson parameter 
A+ = epE(and, analogously, for negative charges). 

Table [3] shows that the collection numbers ccr are very 
small for all considered cases (at least 10 3 times smaller 
than if). Therefore, from Eqs. (fTTil) and (fl2l) we con¬ 
clude that positive charges are not affected by CR col¬ 
lection and have a Poisson distribution for all Z > 2; 
the corresponding Poisson parameter in the region O 
can be as large as A+ ~ 0.4 and ~ 3.1 for the mod¬ 
els 2z? and Jf?, respectively (right panel of Fig. [4j for 
the regions I and C, A+ <C 1). On the contrary, the 
photoemission numbers cpe are usually large (epE f, 
except for the densest C case) and the distribution of 































































































































negative charges is completely dominated by the state 
Z = —1, whose abundance is given by Eq. (fllT) ; for 
\Z\ > 2 it abruptly decreases as Nz oc (\/mecR/epE)^, 
with s/me. cr/cpe ^3x 10 -3 . 

We conclude that in all regions of the cloud, the equi¬ 
librium charge distribution is governed by a competition 
between the cold-plasma collection (providing negative 
charging) and the photoemission due to CR-generated 
UV field (which leads to positive charging). The effect 
of photoemission is dominant in the region O, quite sig¬ 
nificant in the region I (particularly, in the right panel 
of Fig. ©, and still noticeable in the densest region C. 
The direct effect of CRs on dust charging - the collec¬ 
tion of suprathermal CR electrons - is negligible for all 
considered situations. 

3.3. Dependence on the ion mass and grain size 

The results, shown in Fig. © for HCO + ions, remain 
practically unchanged also for other ions: From Eq. © 
it follows that epE oc y/fii, and therefore the charge distri¬ 
bution for Z > 0 [Eqs. (flOl) and (fl2l) ] does not depend on 
rh as long as cpe 1, i.e., for all considered cases. For 
Z < 0 [Eqs. (fTIll and (fl3|) ]. the distribution is unchanged 
when cpe (p. From Table [3] we see that this condition 
is only violated in the densest region C with the “typical” 
parameters (left panel of Fig. ©, where N_i/No oc 

Using the same consideration, one can also obtain the 
dependence of Nz on the grain size: From Eqs. m and 
m we see that, as long as ip 1, the relative abundance 
of the positively charged states varies as Nz/N 0 oc Cp~ x oc 
a. Similar to the dependence on m , the distribution of 
negative charges depends on a only in the region C, where 
epE < <p; in the regions O and X, the size dependence 
sets in when the term e~^ , neglected in Eq. (fl3l) [see last 
Eq. ©], becomes comparable to ccr- The latter occurs 
for sizes > 1.7/(7\o hi £q R ) which are larger than the 
upper cutoff a max of the MRN distribution. 

By employing these simple scaling relations and using 
the examples presented in Fig. 01 one can easily deduce 
the form of Nz for arbitrary plasma composition and 
grain size. 

3.4. Implications 

Knowing the charge distribution on dust grains in 
dense molecular clouds is important for several reasons: 
The charges modify the cross sections of the ion accretion 
on dust, thus critically changing the surface chemistry, 
influencing the formation of grain mantles, etc. For the 
same reason, the charges change the total energy balance 
of grains and, hence, may alter the equilibrium temper¬ 
ature. However, the most pr ofound effect can be on the 
rate of dust coagulation dDominik et al.l [2007 1. Let us 
briefly elaborate on this important point. 

Using the derived charge distributions, one can identify 
the “optimum” balance between the positive and nega¬ 
tive charges, for which the rate of dust coagulation is 
maximized. Conditions for such balance is particularly 
easy to obtain when the distribution of positive charges 
monotonously decreases with Z (i.e., when A+ < 1). 
The coagulation rate in this case is maximized when 
N. |_i ~ iV_ i- For ccr <C 1 and 1 <C <p <C cpe (see Ta¬ 
ble©, from Eqs. (HOD and (fill) we obtain the “optimum” 



Remarkably, the optimum relation becomes size- 
independent at late stages of dust coagulation: When 
the resulting clusters grow well beyond the upper cutoff 
size of the MRN distribution, so that <p < 1, one has to 
replace (p with unity [neglected in Eqs. m and mi 
This yields 

eX/^ ~ E 

Furthermore, since epE oc s/m , the coagulation optimum 
does not depend on the plasma composition either and, 
hence, is universal. As one can see from Table © the 
optimum is attained somewhere between the regions O 
andX, which corresponds to the densities n(H 2 ) between 
~ 10 4 cm -3 and ~ 10 6 cm -3 . 

The existence of the coagulation optimum is in striking 
contrast with the case of pure plasma charging, where 
growing negative c harg es on clusters inhibit their further 
coagulation (lOkuzumil 20091) . Thus, CRs can provide 
ideal conditions for a rapid dust coagulation in cold dense 
molecular clouds. We will discuss this important effect 
in a separate paper. 

4. CONCLUSIONS 

The aim of this paper is to demonstrate that CRs 
strongly affect charging of dust grains in cold dense 
molecular clouds. 

We calculated the local (propagated) CR spectra for 
three characteristic regions of a dense core, and in¬ 
vestigated the two mechanisms of dust charging that 
have been ignored so far: collection of suprathermal 
CR electrons and protons by grains (adding to the cold 
Maxwellian plasma collection), and photoelectric emis¬ 
sion from grains due to the CR-generated UV field. 
While the former mechanism turns out to be always 
negligible, the photoemission is shown to dramatically 
modify the charge distribution for submicron grains in 
the almost entire cloud (as compared to the “reference 
case” of cold-plasma charging). The competition be¬ 
tween the cold-plasma collection (producing, primarily, 
singly charged negative grains) and the photoemission 
(resulting in positive charging) significantly broadens the 
charge distribution. 

The relative magnitude of the CR-induced photoemis¬ 
sion is quantified by the dimensionless number epE- This 
number depends on several physical parameters, some 
of them being only approximately known. The mains 
sources of uncertainty are the values of the ionization 
fraction x e and the CR-induced UV flux X/jv : As we 
pointed out in Sec. 13.21 the uncertainty in x e is mainly 
due to competition between different recombination pro¬ 
cesses, while Xuv, in turn, is determined by (approxi¬ 
mately known) dust albedo, extinction slope, and pho¬ 
toemission yield, as well as by (partially unknown) H 2 
excitation rates. To take all these factors into account, 
in Fig. ©we presented two sets of plots showing the dust 
charge distribution: The left and right panels correspond 
to the expected “typical” and “maximum” values of £pe 
and, thus, demonstrate the extent to which the integral 
effect of uncertainties may affect the final results. 

Our results have several important implications. The 
shown modification of the grain charge distribution con¬ 
siderably changes the rates of ion accretion on dust, 
which in turn can critically change the surface chemistry, 
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alter the total energy balance of grains, influence the for¬ 
mation of icy mantles, etc. These CR effects are partic¬ 
ularly strong in the outer regions of the core, where the 
charge distribution is dominated by positive grains and, 
hence, the accretion of negatively charged ions should be 
drastically reduced. 

The most profound effect of CRs is expected to occur 
for the rate of dust coagulation: When the cold-plasma 
collection is the only charging mechanism operating in a 
cloud, the average (negative) dust charge increases pro¬ 
portionally to the size. Therefore, the growing Coulomb 
repulsion inhibits coagulation of larger (> 1 /im) aggre¬ 
gates. Here we showed that the competition between the 
cold-plasma collection and photoemission can create ap¬ 
proximately equal abundance of positively and negatively 
charged dust, providing “optimum” conditions for coag¬ 
ulations. The derived optimum is size-independent for 
large dust, which enables the growth of big aggregates. 

The presented results are obtained assuming that 
charged dust does not affect the overall charge neutrality 
in a cloud. Although this assumption can be violated for 
sufficiently dense regions, our approach can be straight¬ 


forwardly extended to this case as well, provided the dust 
size distribution is known. Furthermore, the described 
effects of CRs are of generic nature, and are expected to 
operate not only in dense molecular clouds but also in the 
upper layers and the outer parts of protoplanetary discs, 
where mutual sticking of dust aggregates is the essential 
process toward planetesimal formation. 
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APPENDIX 
APPENDIX A 

COULOMB LOSS FUNCTIONS 

In order to combine the propagated CR spectra obtained in lPadov ani et a ll (l2009tl with the cold Maxwellian plasma 
background, the energy loss functions Lk(E) (plotted in Fig. 7 of iPadovani et al.fl2009h have to be extended to lower 
energies (down to about 1 meV), to include Coulomb losses that dominate in this energy range. The Coulomb loss 
term, L^(E), is parameterized for protons bv iSchlickeised (1200211 as 

L°(E) ~ 4.8 x IQ-^XeE- 1 eV cm 2 , 

where x e is the ionization fraction, the proton energy E (in eV) is supposed to exceed ~ 2>k^T. For Coulomb electron 
losses, we used the analytic fit bv lSwartz et all 1 1971 (1. 

L°(E) ~ 5.7 x e E~ 0 - 94 eV cm 2 . 

Dens e cores h ave typ i cal te mperatures of about 10 K, so /cbT - 9x 10 -4 eV. Figure [5] shows the energy loss functions 
from IPadovani et all (12009H with the Coulomb terms included, assuming an average ionization fraction of x e = 10~ 7 . 


APPENDIX B 

SPECIFIC INTENSITY OF SECONDARY ELECTRONS 


After crossing a column density dN = ndx in medium of density n, an electron of energy E and velocity v e has lost 
an energy dE = —L e (E) dN, where L e (E ) is the energy loss function of electrons, 


L e (E) 


1 dE 
n dx 


1 dE 
nv e dt ’ 


(Bl) 


(and similarly for protons). The column density N stop required to stop an electron of the initial energy E, 



dE' 

L e (E'Y 


has a constant value N stop (E) ~ 10 18 cm' 2 up to E ~ 1 keV (see Fig. 8 in lPadovani et al.ll2009f> . and then increases 
weakly with energy. This short range of secondary electrons, in comparison with typical column densities of molecular 
clouds (~ 10 22 —10 23 cm -2 ) justifies a local treatment of ionization. In this “on-the-spot” approximation, the stopping 
time of secondary electrons is 


Tstop (A) 


E 

nv e L e (E ) ’ 


(B2) 


and their stopping range is ~ v e r s t op ~ E/[nL e (E)\. 

Assuming isotropy, the number of first-generation secondary electrons of energy E (per unit energy, volume and 
time) is calculated as 


d/V e sec 
dEdV dt 


= 47m 


h+E 


drr lon 

dE' j k (E')^(E,E'), 
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Fig. 5.— Energy loss functions for protons (black) and electrons (gray). The ionization fraction x e = 10 7 has been assumed to include 
Coulomb losses. 


where I = 15.6 eV is the ionization potential of H 2) j k {E') is the intensity of primary species k, and dcr^/dE is the 
corresponding differential ionization cross section. The number of secondary electrons produced per unit energy and 
volume is then 


dA/| ec 

dEdV 


r sto P = 

dEdVdt v e L e (E) 


h+E 


drr [on 

dE’ 




This quantity is related to the specific intensity of secondary electrons jl ec (E) (number of electrons per unit energy, 
area, time and solid angle) by d/V| ec / dE dV = (47r/ v e )jl ec , which finally yields 

J? r°° drr^ on 

j T{E) * l E dE' j k (E')^(E,E'). (B3) 

Equation (IB3I) is iterated, to compute intensities for the next generations of secondary electrons. 


APPENDIX C 

ELECTRON STICKING PROBABILITY AND SECONDARY EMISSION YIELD 

Following iDraine fc Saloeteil (I1979D . we set the sticking probability s e (E) equal to unity if the electron stopping 
range R e (E) (in the dust material) is smaller than 4a/3, otherwise s e {E ) = 0. The transition energy varies with dust 
size as oc a 2 / 3 ; for sub-micron grains , a ~ 0.1 pm,the transition occurs at E ~ 3 keV. The electron trapping by small 
(PAH) grains was calculated bv lMicelotta et al . (|201 Cf). 

The secondary emission yield S e (E) is obtained by averaging over the velocity distributions of the emitted electrons, 
which are believed to be broad (non-Maxwellian), decaying as E~ x at large energies (IDraine fc Salneterlll979D . For a 
positively charged dust the explicit dependence on ipz is approximated by 8 e oc 1/yT + (< pz/E *) 2 , where F* ~ 3 eV. 
This factor can be neglected using the same reasoni ng as for t he OML c orrection factor omitted in Eq. d8]). For the 
dependence on E we employ the Sternglass formula dHoranvi et aljll988 ). 


Se(E) 

JJmax 

u e 


E 

■S'max 


exp 



where the value of <5™ ax is typically between 1.5 and 2.5, and -E max = 0.2 — 0.4 keV ([Draine fc Salpeterlll979[) . Note 
that S e (E) has a finite-size correction at large E , which increases the yield by a factor of ~ 1.6. 
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